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The concept of incomplete vector-meson dominance and Regge models is applied to the transition 
form factor of the pion. First, we argue that variants of the chiral quark model fulfilling the chiral 
anomaly may violate the Terazawa-West unitarity bounds, as these bounds are based on unverified 
assumptions for the real parts of the amplitudes, precluding a possible presence of polynomial terms. 
A direct consequence is that the transition form factor need not necessarily vanish at large values of 
the photon virtuality. Moreover, in the range of the BaBar experiment, the Terazawa-West bound 
is an order of magnitude above the data, thus is of formal rather than practical interest. Then 
we demonstrate how the experimental data may be properly explained with incomplete vector- 
meson dominance in a simple model with one state, as well as in more sophisticated Regge models. 
Generalizations of the simple Regge model along the lines of Dominguez result in a proper description 
of the data, where one may adjust the parameters in such a way that the Terazawa-West bound 
is satisfied or violated. We also impose the experimental constraint from the Z 7ro7 decay. 
Finally, we point out that the photon momentum asymmetry parameter may noticeably infiuence 
the precision analysis. 
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INTRODUCTION 



The pion transition form factor, F^o^^ , measured in 
the e+e^ — > 7r*'7 annihilation, has been a particularly 
interesting object of study since its value at the origin 
is fixed by the chiral anomaly [H Q while the behav- 
ior at very high Euclidean momenta has been predicted 
via perturbative quantum chromodynamics (pQCD), ap- 
parently holding at sufficiently high but unspecified 
scales [3, 13] ■ However, objections have been raised [1] 
and reiterated @ against the applicability of pQCD to 
exclusive processes. The recent measurement of Fj^o^^, 
by the BaBar Collaboration 0, where the pion tran- 
sition form factor goes visibly above the conventional 
pQCD prediction at scales > 15 GeV^, raised seri- 
ous doubts concerning the theoretical understanding of 
the exclusive processes, also in the asymptotic region. In 
the pQCD approach based on factorization, one uses the 
light-cone Feynman rules and the ERBL evolution of the 
pion distribution amplitude (PDA) 3, 4, 8-16], which re- 
sults asymptotically in the leading-twist Brodsky-Lepage 
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2/^/(5^, in vivid contradiction to the re- 
cent BaBar 0] data. 

Several ideas have been proposed to circumvent the 
problem. Radyushkin [l^ pointed out that the presence 
of the possible end-point singularities in the PDA (as 
found by the authors in the chiral quark models at the 
low-energy quark- model scale (IS]), together with essen- 
tially switched-off evolution and regulated quark propa- 
gators, is capable of reproducing the data in the avail- 
able domain. Similar conclusions were reached by 
Polyakov [l9l ]. In this approach asymptotically i^^o.^^. ~ 
log((5^/At^)/Q^, with the log in the numerator indicat- 
ing the breaking of factorization. We note that the same 
asymptotics follows in the Spectral Quark Model (SQM), 
see Eq. (14.1) of Ref. In fact, serious concerns 

have been spoken out on the validity of factorization [2lj 
in the considered process. Dorokhov [23 - 123 ] proposed 
the use of the fixed-mass constituent quark model to 
evaluate the triangle diagram of Fig. (TJ which is capa- 
ble of reproducing the BaBar data, however the needed 
value of the constituent quark mass is uncomfortably low, 
M 135 MeV. In this model the asymptotics has the 
form Jlog((5^/^^)]^/(3^. Possible need of the higher- 
twist terms has been invoked in Ref. [1^ . The calculation 
of Kotko and Praszalowicz in the non-local chiral 
quark model inspired by the instanton-liquid model of 
QCD produced the result in agreement with the data at 
lower values of Q^. The influence of non-perturbative 
gluonic components of the pion stemming from instan- 
tons has been considered in Ref. [13] . Mikhailov and Ste- 
fanis [l^ showed that the significant experimental growth 
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of the transition form factor between 10 and 40 GeV^ 
cannot be explained in terms of higher-order pQCD cor- 
rections at the NNLO level. Finally, in a very recent 
paper Dorokhov reiterates the findings of Ref. (2^ 
in the non-local model, where the asymptotic behavior is 
- log(QVM^)/Q^- The non-local model of Ref. 28] also 
requires a very low constituent quark mass in order to 
reproduce the data in the whole range. 

Given the fact that the recent BaBar data 0] predate 
the conventional pQCD expectations at such high virtu- 
alities as 15 — 30 GeV^, it seems adequate to question 
any implicit assumptions in theoretical analyses. In this 
paper we analyze critically some of the assumptions un- 
derlying certain high-energy theorems which, if holding, 
would forbid the simple pre-BaBar fits based on the in- 
complete vector-meson dominance (IVMD) [2^ (see also 
Ref. [lO] where this result is upgraded with tiny conse- 
quences for to the muon <? — 2). The theorems in question, 
the Terazawa-West (TW) bounds [3l|, HI], were recalled 
by Dorokhov in Ref. [llj . 

The TW bounds [3l|,|3l] (for a review see, e.g., [Ill) ^'^^ 
central to our discussion; they are reviewed in Sec. IIIII 
Their derivation uses the Schwarz inequality to predict 
that for physical momenta ImF^o^^* — 0{l/yjlp) (the 
first bound). If one further assumes that there are no 
polynomial contributions to the real part of the form 
factor, one may conclude with the help of a dispersion 
relation that Fj^o^^* — 0{1/Q) for space-like momenta 
Q (the second bound). The assumption concerning the 
real part is clearly spelled out in the original works and, 
to our knowledge, has never been credibly questioned. 
We should remind the reader that this was a pre-QCD 
bound, so we may test the validity of the assumption 
by analyzing a particular model where both the anomaly 
is satisfied and the second Terazawa-West bound is vio- 
lated. Actually, a reanalysis of the issue using the quark 
model and corresponding triangle graphs was presented 
in Ref. [sS] as an example case where the expected bound 
behavior holds. Here we will consider another variant of 
the quark model which in fact does not fulfill the assump- 
tions, hence leads to a violation of the second bound. 
The key point is to realize that the anomaly fixes the 
value of the form factor at the origin, while the finite 
Q behavior can be determined regardless of the value at 
(5 = 0. This construction requires a subtraction constant 
which actually implies that at large Q the pion transition 
form factor does not go to zero. These arguments sug- 
gest that pQCD might in fact be computing just the sub- 
tracted form factor. Clearly, the issue cannot be settled 
unless one could go smoothly over all the available en- 
ergy range. Of course, quark models are not QCD itself, 
and it is quite possible that our considerations do not ap- 
ply directly to the QCD analyses. However, the present 
paper unveils a warning that from purely field-theoretic 
reasons an asymptotically non-vanishing pion transition 
form factor cannot be rejected. 

With this finding in mind we analyze in detail the 
predictions of IVMD and the Regge models for Ft^o^^* . 




FIG. 1: The triangle diagram used to evaluate the pion tran- 
sition form factor in quark models, with one real and one 
virtual photon (the crossed diagram not shown). 

Regge models with (infinitely many) tree-level meson and 
glueball exchanges are a realization of the large- TVc limit 
[35| which, unlike the quark models, incorporates con- 
finement. In particular, these models where insightful in 
analyzing condensate issues ]^36 - 44] , form factors [isl . |46| , 
or the nature of scalar states [47|] . 

On purely phenomenological grounds, we will show 
that the data for F^o^^* can be fitted very well in a sim- 
ple vector-meson exchange model with IVMD. In more 
elaborate variants of the Regge approach the data may be 
reproduced (with different choices of the model parame- 
ters) almost equally well with or without the fulfillment 
of the TW bounds. 

There has also been an important interest in the pion 
transition form factor after the proposal of determining 
the rare Z — ^ 7r°7 decay [i^, which probes the time- 
like value — M'^. The experimental bound given by 
the Particle Data Group [4§] implies a large suppression 
as compared to the = anomaly point, but not nec- 
essarily as much as predicted by many authors, where 
T{Z ^ Tr°-f) - 10-"GeV (see, e.g.. Ha]). 

The paper is organized as follows. As a preparatory 
material, in Sec. |TT] we provide some basic kinematics, 
in Sec. IIIII we review the TW bounds, and in Sec. IIVI 
recall the bound from the Z ^ tt^j decay. In Sec. |V] 
we use several variants of the quark model to inquire on 
the validity of the TW bounds. As we will show, these 
models are not capable of describing the data with realis- 
tic model parameters, but pose a warning on the second 
TW bound. We use this insight in Sec. |VT] to propose 
a more realistic Vector Meson Dominance description of 
the data. The analysis is complemented in Sec. IVIII with 
the inclusion of the large Nc motivated infinite tower of 
radial Regge like excitations of the p— meson. The im- 
portant role of photon asymmetry is discussed briefly in 
Sec. I Villi Finally, in Sec. lIXI we list our main conclusions. 

II. THE KINEMATICS 

With the outgoing momenta and polarizations of the 
photons denoted as qi,ei and (72, £2 (see Fig.[T|), one finds 
the amplitude 

r^o^.^.(9i.92) = e;,i.a/3e^e2g"g2-PV7-7*(Q^,^), (1) 

where the pion transition form factor Ft^^*^* depends 
on the total virtuality, Q^, and the photon momentum 
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asymmetry parameter, A, 



<ll-ll 



Q' = -{li + 4). A = -1 < A < 1. (2) 

Equivalently, qf = = -^^^Q^. For large 

virtualities, assuming factorization, one finds the stan- 
dard twist decomposition of the pion transition form fac- 
tor i, 



with 



Nr. 



dx 



(2) 



{x) 



1 - (2.T - 1)2^2 ' 



(3) 



(4) 



involving the pion parton distribution amplitude PDA, 

tf{x). 

Taking the asymptotic form ip{x) — 6a;(l— a;) in Eq. (|4]) 
gives the Brodsky-Lepage asymptotic result Q 



J(2)(^ = l) 



6 



Nef 



(5) 



An analysis of the lowest Gegenbauer moments 02 and 
04 of 1^(2:) has been carried out in Ref. ^51,] and 
fla . [52j | . Higher twists have been analyzed in the frame- 
work of the light-cone sum rules [s^. A direct measure- 
ment of the DA has been presented by the E791 collab- 
oration [H^l- For a concise review on all these develop- 
ments see, e.g., Ref. [S^ and references therein. 

Although theoretically one may take the maximally 
asymmetric case of A = 1, in the experiment this situa- 
tion is not possible to achieve. For instance, in the BaBar 
kinematic setup —q\ < 0.6 GeV^ and — q| > 3 GeV^, sug- 
gesting A '^0.9 — 0.97. As advanced in Ref. [56|, we note 
that the model results presented in the following sections 
become quite sensitive to the precise value of A at large 

, hence in precision analyses the effects of kinematic 
cuts should be incorporated. 

At this point, before undertaking further elaborations, 
an important qualification is in order. Eq. Q provides 
a gauge-invariant high-energy-motivated definition of the 
pion PDA, which formally determines a non-perturbative 
matrix element. A purely low-energy definition of the 
PDA involves link operators between the quark fields at 
different space-time points and depends on the renor- 
malization scale. A non-trivial issue has to do with the 
equivalence of both definitions and the scale at which 
this identification makes sense. We refer to previous 
works [H, [13, for a detailed discussion on this 

delicate and non-trivial issue. An outstanding outcome 
was that in several models 



ip{x) = 1 



(6) 



at the model low-energy scale; the end-point behavior 
would be bended after the QCD evolution to higher scales 



(see Eq. (5.26) of Ref. [5g] for the detailed asymptotic 
form of the PDA near the end-points after evolution). 
While the identification of the PDA through Eq. (jlj in- 
deed requires taking asymptotically large momenta, it 
does not necessarily follow that the asymptotic expres- 
sion itself holds at the so far measured momenta. 



III. TERAZAWA-WEST BOUNDS 

Terazawa 3lj and West [l^ (for a review see, e.g., Q) 
derived unitarity bounds for the pion transition form fac- 
tor. The derivation considers the photon propagator, the 
charge form factor of the pion, and the transition form 
factor. The Schwarz inequality is used to show that 



ImR 



[q^) = 0(l/vV) 



(7) 



for time- like momenta, q^ > 4m2, which we term the 
first TW bound. Under a further important assump- 
tion that there are no polynomial contributions to the 
real part of the form factor, one may also conclude that 
|F^o.^..y.((72)| = 0{l/\/(f). Then the dispersion relation 
(see Sec. 14 of Ref. ) yields that 



\F^o^r'iQ^)\=0{l/Q) 



(8) 



for large space-like momenta Q, which we call the second 
TW bound. 

The assumption of the absence of the polynomial terms 
(in particular, a constant) from the real part of the pion 
transition form factor is equivalent to validity of the un- 
subtracted dispersion relation. As mentioned in the In- 
troduction, this crucial assumption is clearly spelled out 
as such in the original papers; in the following Sections 
of the paper we will investigate situations where the sub- 
traction constant is present in the pion transition form 
factor. 

The constant in the bound ^ may be given [ssl] in 
terms of the photon spectral density and the pion quark 
structure function, namely 



F,o^^.(Q2)| < 



2VnM f\ Fi{x,Q^) 



Q 



where 



His) 



167r3a^ED 



^^V ^(1 ~ ^) 



e — >-hadrons ('^) ; 



with the asymptotic value 



(9) 



(10) 



(11) 



Ci denoting the charge of the quark of fiavor i in the 
natural units. The pion structure function is 
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with Q^) and qi(x, Q"^) denoting the quark and anti- 
quark distribution functions at momentum Q. 

With the help of the SMRS ^ and the GRV 
parameterizations of the pion PDF's we obtain, via in- 
tegrating from X — 10~^ to 1, a numerical estimate for 
the coefficient in the bound © which actually exhibits 
a rather mild dependence on Q^. With in the range 
10-40 GeV^ we find for the LO and NLO evolved PDF's 
the bounds 

|F,.„.(<?=)| < ?S1 (LO), 

< ^ (NLO), (13) 

with the small uncertainty refiecting the considered Q 
interval. Therefore the bound is completely "inefficient" 
in the considered range of momenta, as it is an order of 
magnitude above the BaBar data. For that reason it is 
of formal rather than practical interest. ^ 



IV. BOUND FROM THE Z -> Tv^-y DECAY 

Another bound for the pion transition form factor 
comes from the rare Z — > ir'^j decays. In that process, 
proceeding via a quark loop, only the vector coupling of 
the Z boson to the quarks contributes [i^, hence 



r(0) 



(14) 



The experimental hmit r(Z 7r%) < 5 x 10"5rt,ot(2') = 
10.25 X lO-^GeV, provided by the Particle Data 
Group [4§|, implies 



\Fz^^o^{A4)/Fz^,o^m <0.17 



(15) 



This experimental bound is not as stringent as predicted 
by m any authors, where T{Z — >• 7r°7) ~ 10~^^GeV (see, 
e.g., [501). The bound can be used in models which pre- 
dict the form factor in the time-like region, e.g. the Regge 
models considered in the following Sections. 



V. QUARK MODELS 

The role of the constituent vs current quarks in the 7r° 
and uj transition form factors in triangle diagrams was 
recognized in Ref. [6l| . The diagram of Fig. [T] is superfi- 
cially linearly divergent, but the requirement of the gauge 
invariance guarantees convergence. In the constituent 



^ We are taking large enough as to neglect the higher twist 
corrections. Actually the bound never crosses the pQCD limit 
2/7r/Q^ for the lowest momenta considered in Ref. [60l |. namely 
Qg ~ 0.26 GeV^ at LO and ^ o.40 GeV^ at NLO. 



quark model (CQM) a direct computation of the triangle 
yields, for the general kinematic case. 



1 



(16) 



where the loop function (we work in the strict chiral limit 

of TOtt = 0), 



G{Q\A) 



2AP dx 



(17) 



log 



2M2 + + A)x{l - x)Q2 
2Af2 (1 - A)x{l ~ x)Q^ 



is normalized to unity at the origin, G{0, A) = 1. For the 
special case A — 1 

G{Q^,A = l) = G{Q^)^ 



2M^ 



dx 

— log 



l + x(l-.) — 



(18) 



We note that the evaluation is covariant, thus different 
in philosophy from the conventional light-cone analysis. 
In the latter case the soft part of the diagram containing 
the PDA is factorized, while the hard part of the diagram 
is evaluated according to the light-cone Feynman rules. 
On the other hand, in CQM one evaluates the diagram of 
Fig. [1] with the instant- form Feynman rules, correspond- 
ing to local or non-local variants of the model. That way 
one is insensitive to the issues of factorization. However, 
it is not clear how credible this approach is at high vir- 
tualities, in particular in local models, where the virtual 
quark propagator carries, even asymptotically, a large 
constituent mass. In the study of this Section we are, 
however, primarily interested in formal aspects, namely, 
the TW bounds in a field-theoretic model, thus we shall 
not be concerned with the issue whether or not a low- 
energy quark model can be realistically used to obtain 
the asymptotic transition form factor according to Fig.[T] 

Chiral quark models arc particular realizations of the 
large- A"c limit. A variant of the CQM, the well-known 
Georgi-Manohar model |62|, allows the quarks to carry 
an axial charge different from unity, ^ 1 . The relevant 
part of the Lagrangian of the model is 

L = q(i^ + .g^475 - m) q + ^Tr {^^U^^^'U) + WZW, 
A^^^{u^d^u-ud^u^), u = e**-^/(2/), U^u\ (19) 

where / = 93 MeV is the pion decay constant and WZW 
denotes the Wess-Zumino-Witten term [gsI . \q^ . There 
has been some discussion on whether or not 1 — 5^ = 
0{N^). This issue was answered in the affirmative way 
in Ref. [1^, where the 0{N^) departure from unity is 
basically due to the t-channel exchanges in the Adler- 
Weisberger sum rule for the pion-quark scattering. In 
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the NJL model with vector mesons, the Ai — tt mixing 
is an exphcit leading- A^j, source |66l - [68t of the effect. Be- 
cause of the chiral anomaly, the field representation used 
to describe the 5^ 7^ 1 situation is relevant. In Ref. f69] it 
was shown that the apparent anomalous inequivalence of 
effective theories may be compensated by including extra 
terms in the effective action. More specifically, the origi- 
nal GM model does nof contain the anomalous piece, such 
that the Wess-Zumino-Witten (WZW) action should be 
added a posteriori, as in Eq. (jl9l) . An example on how 
various vertex functions are modified by naively includ- 
ing ^ I and violating the anomaly in chiral quark 
models is presented in Ref. [73|. The subtraction proce- 
dure described above to restore the anomaly and based 
on [6§| was addressed in Ref. [7l| (see also [f^)- 

A direct calculation of the pion transition form factor 
in the GM model supplied with the WZW action yields 
the result 



F 



1 



Q 

9a 



[GiQ') - 1] , (20) 



with G{Q'^) given by Eq. ([T5)) . This expression satis- 
fies the anomaly and the dispersion relation but does not 
vanish at infinity, since 



1 



9 a 



gi^P [\og{QVM^)\ 



+ 



(21) 



fulfills the first TW bound (for the imaginary part) but 
does not fulfill the second bound (for the real part). 
Similarly, within SQM ^20i] a direct implementation of 
is equally possible and one gets the result of Eq. ([20)1 
with 



2ml 



(22) 



With = 1 this model fulfills qualitatively the result of 
Radyushkin with a similar mass scale , since 



-PVo77'(<3^) — 



, 9Wp HiQ^ml)] 



(23) 



In the local quark models we have the general relation 

9i 



F^o^^,{t) 



1-9^ 



dsF^fis) 



[2F^T(0 



(24) 



which correlates the charge and transition form factors. 
This relation shows that in these models even in the case 
g*^ — 1 one cannot have complete VMD simultaneously 
for the two form factors. 



Unfortunately, attempts to fit the experimental data 
with the quark-model formulas, with t/^ equal or different 
from unity, are not numerically successful in the whole 
momentum range, unless one uses unrealistic model pa- 
rameters. The purpose of the above calculations was dif- 
ferent. Our examples show in an explicit manner that the 
second TW may be violated in models consistent with 
all field-theoretic constraints, such as covariance, gauge 
invariance, chiral symmetry, and anomaly matching. Im- 
portantly, they provide instances where the form factor 
does not vanish asymptotically. The calculation also il- 
lustrates that the — value and the finite values 
are independent. The possibility of disobeying the second 
TW bound will be explored in the following Sections. 



VI. VECTOR MESON DOMINANCE MODELS 

We now come to the core of our paper and consider 
several implementations of I VMD, both in a simple ap- 
proach with a single vector meson, as well as in more so- 
phisticated Regge models with infinitely-many radially- 
excited states. It will turn out that the pion transi- 
tion form factor can be accurately described with this 
phenomenological method. The coupling of photons to 
vector mesons has a long history and Ref. [7^ compre- 
hensively reviews the interplay between universality, the 
vector-meson dominance and the low-energy theorems. 
We will highlight first the issue of IVMD for the charge 
form factor, such that our points are later on more easily 
made for the transition form factor. 



A. Charge form factor 

The charge form factor is the famous case where VMD 
can be implemented: 



Km 



M 2 - t 



Ml 
t 



(25) 



with t = -Q^. When the unregularized quark-loop mean 
squared radius [74;] is matched to ([^5)) . one gets the rela- 
tion 



M2 = 2An^f^/N, 



(26) 



There is no way to match the successful form factor of 
Eq. (UHl) to the pQCD result.^ 



^ If we nevertheless match Eq. I|25|l to pQCD, the only possible 
solution is as{Q) = 't/2, which yields a too small scale Q ~ 
300 MeV. This matching to pQCD may seem weird but need 
not necessarily be conceptually wrong. Note that in the string 
model calculation of the qq potential one gets Vqq{r) = — 7r/12r, 
whereas from pQCD Vqq{r) = —4as/3r. This yields a = 7r/16, 
which means a safely high scale of ^ = 2 GeV for Agf^jj = 
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If one uses a once-subtracted dispersion relation and 
imposes the current conservation, F°"^{0) — 1, one gets 



which yields 



— t — ie 



dt', 



(27) 



F-," (-0O) -^--If ^"^if^^'^ dt'. (28) 
Saturation with a single resonance gives the result 



2 M 2 - t 



a 
2 



a Af2 
2~t 



(29) 



When the coupling a/2 = fp-ygpTvw/ml = 9 p-ktt 1 9 pe+ e~ 
becomes unity, the form factor vanishes at infinity. Saku- 
rai's universality indeed requires (7p,r7r = gpe+e- ■ In gen- 
eral, however, we may have gp,r7r 7^ gpe+e- ■ Note that 
positivity requires a < 2, such that F^™{—Q^) < 1. A fit 
to the experimental data [75l - l7i | yields a/2 = 0.99(2) and 
My = 679(36) MeV. Thus the data for the pion charge 
form factor are consistent with the complete VMD, but 
certainly do not preclude IVMD, with a/2 departing 
slightly from unity. 

Caldi and Pagels [t^ obtained a similar expression as 
Eq. (211) for the pion form factor from a direct photon 
contribution and a momentum dependent p — 7 vertex 
(see also Ref. [lOl). The discussed properties are nicely 
displayed by the hidden-symmetry approach by Bando et 



8l| , where there is a contact piece and the vector me- 



son term, which dominates completely when the KSFR 
relation is fulfilled. Actually, in Ref. [83 the equivalence 
of this approach to more conventional ones is established. 
The interplay between VMD and universality was ana- 
lyzed by Schechter ^SSj] (see also Ref. [13 )■ Symmetry 
breaking effects have been analyzed in Ref. [83] • More 
attempts including predictions for meson decays can be 
found in [86| . 

Of course, the non vanishing of the low-energy rep- 
resentation (the one-resonance saturation) need not be 
taken as a fundamental problem. The only feature we 
see is that this non- vanishing represents more accurately 
the unknown high-energy data. If we had infinitely many 
states, we could fit that data and by separating explicitly 
the contribution from the lowest p-state. We see that the 
effect of all other states does behave as the a constant, 
which is slowly dependent on Q in a wide energy range. 



Ideally, one should take the e+e 



TT+TT data over 



all possible momenta. However, obviously the experi- 
mental data are available only up to a certain maximum 



240 MeV. The string model describes accurately the lattice data 
in the short distance region. In the case of the charge form factor 
with VMD, the required scale is much smaller. 



value, Smax = 4A2. Thus, even if the form factor vanishes 
at infinity, we have 



TT J4„2 S S-t 

t luiF^fis) 



1 



TT 74A2 S S ~ t 



ds, (30) 



where the last term is weakly momentum-dependent and 
hence resembles a constant behavior assumed in IVMD, 
as discussed above. 



B. Transition form factor 

For the tt'^jj* form factor the complete VMD with just 
one state implies 



Ml 



47r2/ Af2 - t 



Ml 



(31) 



Note that while the anomaly value does not depend on 
Nc explicitly, the high momentum behavior does. When 
matching Eq. (I5T]) to the pQCD result of Eq. (O is done, 
one gets independently the relation p6|. 

Ml = 247rV'/^c. (32) 

Despite this appealing property, the parametrization (pij) 
fails to describe the experimental data in the high mo- 
mentum region. 

If we incorporate the possibility that the form fac- 
tor need not vanish at infinity, we may write a once- 
subtracted dispersion relation [87[, 



-(0) 



1 r tlmF^^it') 



t' t' ~t-ie 



dt'. (33) 



The infiuence of the well-known time-like region does not 
determine unambiguously when the onset of the pQCD 
takes place. Actually, the single VMD model shows that 
even in the space-like region with momenta as low as 
~ the effects of the chiral logs and final state 
interactions are meager. Taking the limit t — ^ — 00 we 
get 



to 



t' 



This shows that 



F^Q^^.(-cx)) 



< 



(34) 
if 



Im_FVo^^. (t) > 0. When we saturate the absorptive part 



jusi 

F^o 



with just one resonance, we get 

1 



7r^77 



4^2/, 



1 



'Ml 



(35) 



with My denoting the vector-meson mass. The coeffi- 
cient c is related to the p — decay. 



-2efp Ajp" ^ ^°7) 
nip A{Tr'^ 77) 



1.022 ±0.051, (36) 
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FIG. 2; The pion transition form factor in the space-like re- 
gion = —t > 0. The dots, squares, and diamonds corre- 
spond to the CELLO CLEO Ull, and BaBar Q data, 
correspondingly. The line is the incomplete-vector-meson- 
dominance fit with formula (1351) and parameters (|38p . 



where the value of the estimate is obtained from the de- 
cay width T{p+ ^ 7r+7) = 68 ± 7KeV. 

If we fit the parameters in formula (|35t with the CLEO 
data only, we get 



c = 0.998(18), My = 777(44) MeV, 



(37) 



with x^/DOF = 0.54, hence c is consistent with unity, in 
agreement with previous determinations assuming com- 
plete VMD, however not excluding IVMD. On the other 
hand, the fit to the combined CELLO, CLEO, and BaBar 
data yields the result 



c = 0.986(2), Mv ^ 748(14)MeV, 



(38) 



hence c is different form unity at the level of 7 stan- 
dard deviations, while it remains consistent with estimate 
([55)1 . As seen from Fig. [21 the agreement with the experi- 
mental data over all the momentum range is remarkable, 
with x^/DOF = 0.7. The contours of the confidence- 
levels are displayed in Fig. [3] 

The logarithmic slope at the origin reads 



1 



r-F'7r«77*(Q) 



(39) 



Numerically, taking the optimum parameters psp , we get 
the value 



^ = 1.76(7) GeV- 



(40) 



Our IVMD model estimate is in a very good agree- 
ment with the average experimental values quoted in the 



PDG 



= (1.76 ± 0.22)GeV" . This value is very 



close to that originally reported by the CELLO collabo- 
ration [1^, obtained from an extrapolation from high-Q^ 
data to low by means of generalized vector meson 
dominance, = (1.4 ± 1.3 ± 2.6)GeV"^ given in [93, 
and = (1.4 ± 0.8 ± 1.4)GeV-^ given in 




FIG. 3: The Ax^ = 2.3 and 4.6 contours, corresponding to 
the 68% and 90% confidence levels, in the Mv ~ c plane for 
the fit with the IVMD ansatz ((SSj. The central values yield 

xVdof = 0.6. 



To end this Section, we turn to the rare decay Z — > 
7r'^7, which probes the transition form factor in the phys- 
ical region. From Eq. ([55)1 with parameters psp we get 
\Fz^.^o^{Ml)/Fz-^^o^{0)\ = 0.014(2), a comfortable or- 
der of magnitude smaller than the upper experimental 
bound (dni) , but not as small as predicted by several mod- 
els (~ 10~^, see, e.g., (s^l and references therein). 



VII. REGGE MODELS 

In the previous Section we have considered the sim- 
plest possible implementation of IVMD, with just one 
vector-meson state. However, the large- Vc limit of QCD 
involves tree-level diagrams with infinitely many states, 
including the radial excitations. That way the matching 
to the QCD correlators can be accomplished [s^, H^, |43| , 
as well as the correct asymptotic behavior of the pion 
charge form factor [i^ may be obtained. In this Sec- 
tion we analyze the pion transition form factors in the 
framework of Regge models with infinitely many radially 
excited vector-meson states. 

Based on the success of the Veneziano-Lovelace- 
Shapiro dual resonance model (see, e.g., ^95, '93'] and ref- 
erences therein) Suura [94] and Frampton [95] proposed 
analytic models. The general form reads 



E 



(9? 
92), 



Ml 



(41) 



where Fy^ and Fy^ are the current-vector meson cou- 
plings, while Gt^VpVu, is the coupling of two vector mesons 
to the pion. The situation is depicted in Fig. U) At the 
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v 



FIG. 4: The pion transition form factor in Regge models. The 
labels Vp and denote the infinite vector meson towers with 
the p and uj quantum numbers. 



soft photon point, corresponding to the neutral pion de- 
cay 7r° — 27, the chiral anomaly matching condition im- 



poses the normalization 



VpVuj 

1 



2Fy^(0)FyJO)G,V^yJO,0) 



47r2/ 



(42) 



This consistency constraint, realized in nature, can be al- 
ways satisfied in models by an appropriate choice of the 
couplings. In Ref. 45] we analyzed the transition form 
factor with the help of formula (|4ip with the coupling 
constants Fy taken as constants (as requested asymp- 
totically by the matching of spectral densities to QCD) 
and allowing for constant diagonal couplings in G^VpV^. 
As a result, a constant PDA of Eq. [S]was extracted. The 
corresponding pion transition form factor has the asymp- 
totic behavior ^ \og{Q^ / ji^) / Q"^ , however with physical 
values of the model parameters it overshoots the data in 
the region above ^ 1 GeV^, thus does not properly 
reproduce the data.'^ 

In the analysis below we use a more sophisticated ap- 
proach, namely the factorizable product suggested by 
Dominguez [96j in the dual description of radiative de- 
cays. 



F^O^^* (^) 



1 



Mt), 



where 



(43) 



(44) 



E 



r{2 -b + n) 



r{n + l)r(2 -b)an + - t' 



The function fb{i) depends on three parameters: the 
lowest-lying vector-meson mass. My, the string tension. 



^ As already mentioned, we stress that taking the — > 00 limit 
is an operational way of extracting the leading twist PDA; this 
is different than describing the current data within this limit. 



(7 = a/(27r), and the asymptotic fall-off parameter, b. 



function (l45l) fulfills the normalization condition 



The mass formula is then Mv{nf 



MO) - 1. 



The 



(45) 



For X ^ y one has B{x,y) ^ r(y)a; y, hence in the 
asymptotic region of My — <^(6 — l)awe find 



hit) 



^ + b-i 



6-1 



(46) 



The TW bounds are satisfied if 6 > 1.5. 

We remark that in Ref. [i^ this version of the Regge 
approach was used to describe the charge form factor. We 
have shown that it can be accurately reproduced up to 
~ 6GeV^, while the pQCD result greatly undershoots 
the experiment. The onset of pQCD occurs at extremely 
high ( "cosmological" ) values of Q^. 

We now proceed with the application of the presented 
Regge model to the pion transition form factor. Taking 
a = 1.3 GeV^ (which correspond to the string tension 
a = (455 MeV)2), a fit of formula (gS]) to the joint 
CELLO, CLEO, and BaBar data yields 



My = 0.672(25) GeV, 6=1.81(3), 



(47) 



with x^/DOF = 1. The fit is shown in Fig. [5] with 
the dashed line. Since b > 1.5, the fit satisfies the 
TW bounds. For the logarithmic slope at the origin 
we find b^r = 2.1(2) GeV~^. We also get, with smooth- 
ing the spectral density in the physical region, the ratio 
\Fz-,^o^{M^)/Fz^^o^{0)\ = 0.0014(4), which is com- 
fortably below the experimental bound and an order of 
magnitude smaller than the IVMD fit of Sect. lyTBl 

We can further try to improve the agreement with the 
BaBar data by explicitly separating the first pole as fol- 
lows: 



F 



,{t) = 



1 



47r2/. 



t 



, hit) 

m 



, (48) 



where fbit) is obtained from Eq. (j^H) with the n = 
term omitted from the sum. We impose the TW bound 
by setting b — 1.5. Then (with a = 1.3 GeV^) we find 

c = 0.95(2), My = 709(9) MeV^ (49) 

with xVDOF = 0.85, b^ = 1.9(1) GeV"^ and 
\Fz-,^o^{M^)/Fz^^o^iO)\ = 0.01(1). The correspond- 
ing curve is presented in Fig. [S]with the dotted line. 

Finally, we take a subtracted Regge model, in analogy 
to the model of Sec. IVIBj of the form 



'it) 



1 



47r2/. 



[l^c + cfbit)] . 



(50) 



Setting Mv = 770 MeV and a = 1.3 GeV^ yields the 
optimum values for the remaining parameters 



c = 0.984(4), & = 2.05(3), 



(51) 
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FIG. 5: The pion transition form factor in the space-like 
region = —t > 0. The dots, squares, and diamonds 
correspond to the CELLO [H, CLEO [H, and BaBar 0] 
data, correspondingly. The dashed line is the results of the 
Regge fit with formula (|43[) and parameters ((47]). The dotted 
line shows the Regge fit with Eq. (|48|I and parameters (|49|l . 
Finally, the solid line corresponds to the subtracted Regge 
model of Eq. (|50p and parameters (|51|) . 
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FIG. 6: Same as Fig. [5]for a smaller range of . 
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FIG. 7: The IMVD fit of Eq. ^ with A = 1, 0.975, and 0.95, 
denoted with solid, dashed, and dotted lines, respectively. 



VIII. ASYMMETRY PARAMETER 

As mention in the Introduction, in the BaBar kine- 
matic setup —q\ < 0.6 GeV^ and — > 3 GeV^, sug- 
gesting A ~ 0.9 — 0.97, hence A is not strictly 1. This 
departure has significance for the fits and the obtained 
parameters, which should not be forgotten in precision 
analyses. We take as an example the IMVD model, which 
now becomes 

F^o^^,{-Q^) = X (52) 



Fitting with A — 1, 0.975, and 0.95 yields, respectively, 
c 0.986, 0.978, 0.974, and My = 748, 754, 768 MeV. 
We note significantly diff'erent values of the optimum 
parameters, with My closest to the physical value for 
A = 0.95. The corresponding curves are displayed in 

Fig.m 



with x^/DOF = 0.7. Since the value of b is consis- 
tent with 2, the model gives very similar results to the 
model with the single vector-meson state considered in 
Sec. lyTBl Hence we find = 1.69(1) GeV"^ and 
]Fz^^o^{M^)/Fz-,^o^{0)\ = 0.021(4). The correspond- 
ing curve is presented in Fig.[5]with the solid line. Within 
the present model we see that the higher radially excited 
states of the vector mesons, p', p", . . . and w', w", . ■ • , are 
weakly coupled. 

We also zoom the low-Q^ range in Fig. |6l All the 
considered fits practically overlap in the displayed range, 
in agreement with the observation that IVMD could not 
be rejected even at ~ 8 GeV^. Our analysis shows 
that higher energy data do in fact favor IVMD at the 
level of 7 standard deviations. 



IX. CONCLUSIONS 

These are our main findings: 

• The Terazawa-West bound, asserted for the real 
part of the pion transition form factor, need not be 
fulfilled in field-theoretic approaches. We provide 
an explicit counter-example with the low-energy 
Georgi-Manohar model, where the chiral anomaly 
is fulfilled but the transition form factor does not 
vanish at t —oo. Provided this feature holds in 
QCD, it opens a possibility of explaining the BaBar 
experimental data with models incorporating the 
incomplete vector-meson dominance. 

• Moreover, the coefficient in the Terazawa-West 
bound, estimated with the help of a phenomenolog- 
ical parametrization of the pion parton distribution 
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functions, is large, extending an order of magnitude 
above the BaBar data. Thus, even if it holds in 
the absence of polynomial contributions to the real 
part, it is completely ineffective for the momentum 
range of interest. 

• Already the simplest model with the incomplete 
vector-meson dominance, incorporating a single 
vector-meson state, is capable of reproducing the 
data in the whole available experimental range, 
< Q2 < 35 GeV^ 

• Within the Regge approach, where infinitely many 
radially excited states are included, the data can 
be fitted both ways: satisfying or violating the TW 
bound. The agreement with the experiment is sat- 
isfactory, both near Q"^ — 0, where the anomaly 
value and the slope of the form factor are repro- 
duced, as well as in the intermediate CLEO range 



and the high-Q^ BaBar range. 

• An additional constraint on the models in the time- 
like region of momenta follows from the rare Z — >■ 
7r°7 decay. We use this bound in our considera- 
tions. For the considered models it is comfortably 
satisfied. 

• Finally, we note that the numerical fits are quite 
sensitive to the photon momentum asymmetry pa- 
rameter, A, which leads to sensitivity in the phys- 
ical parameters, such as the vector meson mass, 
and sensitivity of the transition form factor in the 
asymptotic range. Since A is not strictly 1, the ef- 
fect of the kinematic cuts should be considered in 
precision analyses. 
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